Design strategies for the creation of aperiodic nonchaotic attractors 
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Parametric modulation in nonlinear dynamical systems can give rise to attractors on which the 
dynamics is aperiodic and nonchaotic, namely with largest Lyapunov exponent being nonpositive. 
We describe a procedure for creating such attractors by using random modulation or pseudo-random 
binary sequences with arbitrarily long recurrence times. As a consequence the attractors are geo- 
metrically fractal and the motion is aperiodic on experimentally accessible timescales. A practical 
realization of such attractors is demonstrated in an experiment using electronic circuits. 
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Generating dynamics which is aperiodic, but 
nevertheless stable in a sense that nearby tra- 
jectories coalesce and synchronize, has been of 
considerable interest in the past few years. Such 
motion has been typically observed in driven dy- 
namical system, in particular when the drive is 
quasiperiodic. Since quasiperiodicity is difficult 
to achieve in practice, a major issue in this re- 
gard has been whether such dynamics can be 
achieved by other techniques, and has been an- 
swered with various degrees of success. Here we 
present a simple design scheme that uses pseudo- 
random binary sequences with very long recur- 
rence times to switch the dynamics between two 
different states. The resultant dynamics goes to 
an attractor which is aperiodic and stable, namely 
has negative Lyapunov exponent. Characteriza- 
tion of such dynamics reveals the fractal nature 
of such dynamics and also their differences with 
the ones obtained by quasiperiodic drive. Such 
a design scheme is further realized in an experi- 
mental setup using electronic circuits, suggesting 
potential applications in practical situations. 



I. INTRODUCTION 

The design of dynamical systems in which the motion is 
both aperiodic and stable has been an objective of a num- 
ber of recent studies [1]. One realization of this goal is 
in strange nonchaotic attractors (SNAs) that can be cre- 
ated in quasi-periodically driven systems [3, 0] . In such 
systems which have been known for over twenty years 
now [3], the attractor has a fractal geometry, and this 
results in the dynamics being aperiodic. The Lyapunov 
exponents of the drive system are nonpositive: this lead 
to a lack of sensitivity to initial conditions, and thus to 
the synchronization of orbits [4(. 

Achieving quasiperiodicity is simple in principle, but 



difficult in practice. Quasiperiodic modulation requires 
that the system be driven either through a single source 
that has an irrational frequency (for maps), or with 
two sources whose frequencies are incommensurate (for 
flows). Experimental uncertainties are usually larger 
than the precision with which rational or irrational num- 
bers can be measured, and therefore finding or creat- 
ing SNAs in practical situations has proved to be diffi- 
cult, except for isolated experiments [5|]. If strange non- 
chaotic dynamics is to be taken seriously, it is necessary 
to ask if such behaviour can arise in the absence of strict 
quasiperiodicity. In particular, most physical, biologi- 
cal and engineering systems will not be quasiperiodic, so 
it is natural to ask whether SNAs can arise in situations 
where the underlying system is autonomous, or has other 
time dependence. 

Earlier studies have addressed this issue @, 0, H, 
with varying degrees of success. The use of an external 
noise source typically has the effect of smearing out at- 
tractor structure and gives effective Lyapunov exponents 
that are nonnegative. As a result the dynamics is neither 
truly strange nor truly nonchaotic [TO]. The recent sug- 
gestion by Wang et al. [9], that additive noise alone can 
be used to induce robust SNAs in both maps and flows 
appears to ensure that the attractors so created share 
the mathematical properties of SNAs formed by other 
bifurcation routes [l|, 0] • 

The approach taken in the present paper pursues a dif- 
ferent route. We show that by using parametric modula- 
tion based on deterministic pseudorandom dynamics, it 
is possible to create dynamical attractors that are stable, 
nonchaotic and aperiodic, and which will appear strange 
on any measurable timescale. In this context the term 
'strangeness' refers purely to the geometry, while stabil- 
ity refers to the fidelity in signal reproduction or tracking 
[ill ], an issue that, for instance, underlies schemes for en- 
cryption and communication. 

Our procedure for the creation of such aperiodic non- 
chaotic attractors (AN As) relies on the use of binary ran- 
dom numbers to modulate the dynamical system between 
two dynamical states, a stable fixed point and a chaotic 
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attractor. By suitably designing the drive dynamics, it 
is possible to ensure that the asymptotic dynamics of 
the driven system has (a) a negative largest Lyapunov 
exponent, and (b) nontrivial and complicated geometry 
on spatial scales that are determined by the (essentially) 
experimental resolution. We show that such attractors 
can be created in both discrete-time autonomous maps as 
well as flows, and further present an experiment based on 
electronic circuits to support our findings. This suggests 
that ANAs could have potential application in practical 
situations where aperiodic dynamics is desirable, as for 
example in chaotic communications. At the same time, 
these attractors have differences from those created by 
quasiperiodic driving and other methods [H, 0] . It should 
be noted that dichotomous driving has been used before 
in both numerical [l2[ as well as experimental [l3[ stud- 
ies, although the motivation there was to study noise- 
induced transitions between different states or attractors. 

In the following section we discuss the design of ANAs 
in the driven Henon map and the driven Lorenz system 
using a deterministic feedback shift register to generate 
a pseudorandom drive [l| . In Sec [Til] we discuss similar 
driving mechanism using chaotic sequence from a Chua 
circuit. The study of such attractors and their proper 
characterization is discussed in Sec IIV1 where we show 
the geometric differences between ANAs and comparable 
SNAs. In Sec [V] we present an experimental realization 
of such attractors using electronic circuits, and conclude 
in Sec I VII with a discussion and summary of our results. 



II. DICHOTOMOUS DETERMINISTIC 
MODULATION 

Consider a dynamical system (with one freedom and 
a single parameter for simplicity, but with obvious ex- 
tension to higher dimensions and to the case of several 
parameters) 

s -/(*,&), (1) 

that is modulated through the output of a binary drive 
sequence (strings of O's and l's) z n as 

Xn+i = f(x n , h + z n (b 2 - h)). (2) 

Depending on the value of z n , the system parameters thus 
switch between b\ and b 2l giving a dichotomous modula- 
tion that is, nevertheless, deterministic. 

One standard way of achieving this is to use a lin- 
ear feedback shift register (LFSR)pJ| that generates a 
pseudo-random bit sequence {z} through a delay map- 
ping of the general form 

N 

Zn+i = ^CLiZn+i-i mod 2, (3) 



where a is also a binary variable. For a specific choice 
of nonzero a's for a given N (the "tap sequence"), the 
dynamics is on an attractor with period < 2^-1. The 
analog generalization (namely, the analog feedback shift 
registers (AFSR)) [15] uses the same coefficients in a con- 
tinuous mapping 

1 1 N 

Z n +l = - - -COS7T^ai2 n+ i_i, (4) 
i=l 

to generate a pseudorandom sequence of O's and l's. This 
dichotomous drive is, importantly, a dynamical system 
and the drive sequence is an attractor of the dynamics. 
The sequence is optimal if arbitrarily long sequences of 
either or 1 occur. The theory of LFSRs (and thus 
of AFSRs) is well-developed and minimal tap sequences 
that produce the longest possible (namely 2^-1) period 
pseudorandom sequences are easily available [14]. For 
sufficiently large AT, the period of the pseudorandom se- 
quence can quickly exceed the age of the universe at any 
realistic sampling rate. 

Designing aperiodic but nonchaotic dynamics in x is 
straightforward: for instance if b\ corresponds to a case 
of, say, superstable dynamics, and b 2 to the case of 
chaotic dynamics in the system (Eq. [2j), the resultant 
dynamics in the driven system will be aperiodic but will 
rapidly be attracted to the superstable orbit whenever 
there is a "gap", namely a string of O's. The Lyapunov 
exponent will consequently be negative and as a result 
trajectories with arbitrary initial conditions will synchro- 
nize. 

Special attention should be given when choosing the 
tap sequence TV. It must be large enough such that the 
recurrence time, namely 2^ — 1 is much longer than the 
time-scales used for simulations. For smaller N the re- 
currence also become short, and since the AFSR dynam- 
ics is periodic, the system dynamics will go to a periodic 
attractor. 

Any other random sequence will also serve the purpose, 
but LFSRs or AFSRs offer a practical advantage over 
other pseudo random number generators (PRNGs). The 
shift registers are maximally stable [15]: being attractors 
of the dynamics their stability and controllability — unlike 
that of stochastic sequences or PRNGs — is more easily 
ensured. Furthermore, since feedback shift registers are 
dynamical systems as well, the entire drive-response sys- 
tem can be represented as a delay dynamical system. 

We discuss representative examples below. 



A. Henon Map 

The Henon map [l6[ is a well studied two dimensional 
iterative dynamical system, given by the following equa- 
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FIG. 1: (a) Variation of the Lyapunov exponent with a at 
(3 — 0.3 for the unmodulated Henon map. (b) Aperiodic 
nonchaotic orbit in the Henon map, parametrically modu- 
lated through an AFSR of order N=24. The attractors are 
obtained for parameter switching between a = —0.14 and 
a = —1.2. The largest Lyapunov exponent is Ai = —0.16. 
The recurrence time for the AFSR is 2 24 — 1. 
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FIG. 2: Aperiodic nonchaotic attractor obtained in the 
phase-space for the Lorenz system under the AFSR modu- 
lation. Dynamics obtained for parameter values p — 7 and 
c — such that p switches between the values 7 and 30. 
Switch duration time is chosen to be r = 10, and the largest 
Lyapunov exponent is Ai = —0.00985. 



Xfl-^-~\^ l | GiX y t 

Vn+1 = 0X n . 



(5) 
(6) 



For certain choice of (a, 0) and for a certain range of ini- 
tial condition, the system exhibits chaotic or fixed point 
or periodic behaviour (see Fig. Hk). 

Applying the strategy discussed above produces attrac- 
tors which are both nonchaotic and essentially strange. 
In the Henon map, at a = —1.2 the dynamics is chaotic 
and at a = —0.14 the dynamics goes to a fixed point (see 
Fig. Hk). Combining Eq. (jH) and Eqs. (06|), we get the 
following equations for the drive-response system 



Vn+1 



1 + a(l + cz n )xl +y n: 
Px n , 

i N 

-[1 - cos(7r^aiZ„ + i_ 



(7) 
(8) 

(9) 



1=1 



^r, then de- 



If we choose a = —0.14, /3 = 0.3 and c 
pending on whether z n is or 1, the quantity a(l + cz n ) 
will take values of either —0.14 or —1.2. As a result the 
dynamics hops between two different states such that the 
global dynamics is stable and nonchaotic. The Lyapunov 
exponent for the global dynamics is given by A = —0.16. 
Fig. [U(b) shows the attractor in phase space. Clearly 
the attractor looks geometrically strange; the dynamics 
is aperiodic and nonchaotic. 



a modulated Lorenz system [171 ] . 

x = a(y-x), 
y = (l + cC(t))px 
i = xy — /3z. 



xz, 



(10) 

(11) 

(12) 



which has the parameter p changing in a time-dependent 
manner through the variable £, 



tit > t > (n — 1)t. 



(13) 



As in the mapping in Sec III AI above, z n is the output of 
an AFSR (Eq. [4j) and the switch duration r is an addi- 
tional parameter in the problem: it is the time for which 
a trajectory is switched into either of the states. In the 
present problem we have taken r to switch into either of 
the states to be equal, but one can choose mismatched 
t's. 

The attractor shown in Fig. [2] is the result of the dy- 
namics hopping between parameter values p = 7 and 
p = 30. At the latter value the attractor is chaotic 
with the characteristic butterfly structure about two 
symmetric unstable fixed points {—8.79,-8.79,29} and 
{8.79,8.79,29}. For p = 7 the system has two symmet- 
rical attractive fixed points at {—4,-4,6} and {4,4,6}. 
When the parameter switches between the values, the 
dynamics alternates between the stable and the unstable 
fixed point dynamics, resulting in the structure visible in 
Fig. El 



B. Lorenz system 

The same strategy can be applied to a flow: switching 
the dynamics between a fixed-point or limit-cycle and a 
chaotic attractor can result in such strange dynamics in 



III. CHAOTIC MODULATION 

A sequence generated from a chaotic signal can also 
generate aperiodic nonchaotic dynamics. For instance, 
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FIG. 3: (a) The Chua double-scroll attractor in the x — z 
plane. The fixed parameter values are ci =9, C3 =0, mo = 
— f , mi = |. The double-scroll is obtained by tuning the 
parameter C2. Here we take c<i — 14.141. (b) Corresponding 
bit sequence obtained assigning all points with positive value 
of x as '1' and with negative values as '0'. 

consider the Chua oscillator [18[ which is described by 
the following sets of equations 

x = c x {y - x - g{x)), (14) 
y = x-y + z, (15) 
z = -c 2 y-c 3 z, (16) 

where g(x) is given by, 

mix + mi — mo if x < — 1, (17) 
g(x) = rriQX if — 1 < x < 1, (18) 

mix + mo — mi if 1 < x. (19) 

For certain choice of the parameter C3, a chaotic "dou- 
ble scroll" attractor is obtained (see Fig.[3fci). We extract 
an indicator sequence of binary numbers from this chaotic 
attractor (see Fig. 03d) using the prescription that when- 
ever the trajectory is on the left (resp. right) scroll, the 
indicator sequence is taken as 1 (resp. 0). 

This gives a drive signal which is piecewise constant, 
which upon application to the Lorenz system discussed 
in Sec III Bl yields an attractor (see Fig. [4j) that is very 
similar to that obtained via deterministic dichotomous 
driving. The main feature that both these drive signals 
share is that they have long periods when the drive is on 
the stable attractor, and this suffices to ensure that the 
eventual dynamics is aperiodic, and that the attractor 
has a complicated geometry and a nonpositive largest 
Lyapunov exponent. 




X 



FIG. 4: Aperiodic nonchaotic attractors obtained in the 
phase- space for the Lorenz system using chaotic bit sequence 
from the Chua double scroll. The parameter values are p — 7 
and c = such that p switches between the values 7 and 28. 
Switch duration time r is chosen to be unity. The largest LE 
is Ai = —0.03. The other parameters a and f3 have values 10 
and I respectively. 



In the following section, we characterize these attrac- 
tors via measures used in the study of SNAs. 



IV. CHARACTERIZATION OF ANAS 

Computation of the largest Lyapunov exponent shows 
that these attractors are nonchaotic. However, since the 
attractors are geometrically strange only over finite reso- 
lution, they differ in the qualitative and quantitative as- 
pects of their local fluctuation properties [19[ from other 
similar attractors such as SNAs. 

Finite-time Lyapunov exponents (FTLES) [20[ are lo- 
cal estimates for the rate of divergence between nearby 
trajectories, and explicitly depend both on the time inter- 
val r over which they are measured as well as the initial 
conditions. By computing A r for a large number of initial 
points in the phase space, one can obtain the stationary 
distribution, 

P(A, r) = Probability that A r (20) 
lies in the interval (A, A + dX). (21) 

For SNAs this is typically broad and non- Gaussian [H, 
I21L I22I ]. although the mean of the distribution, namely 
the asymptotic Lyapunov exponent, is negative. 

The FTLE distribution for ANAs is purely Gaussian: 
see Figs. 13(a) and (b)) for the driven Henon and Lorenz 
systems. This follows from the pseudo-random nature 
of the driving which switches the dynamics between a 
regular and a chaotic state. Correlations die out rapidly 
and the FTLEs satisfy the central limit theorem. As a 
result the distribution is a Gaussian whose spread is a 
function of the length of the trajectory [19[. 
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FIG. 5: Finite time Lyapunov exponent (FTLE)s. (a) for the 
modulated Henon map with parameter values a = —0.14 and 
c = Asymptotic largest LE value is Ai = —0.16. The 
distribution is taken for N = 100. (b) Similar results for the 
modulated Lorenz system with parameter values as in Fig [2] 
and asymptotic largest LE is Ai = —0.00985. 



The contrast between ANAs and SNAs is also evident 
in the Fourier spectrum Q . Upon computing the time- 
dependent partial Fourier sum [2, [H| 

N 

X(Q, N) = ^2x n exp i27rnr \ (22) 

n=l 

at irrational frequency Q the golden mean ratio, — 
l)/2, the scaling relation |X(^,7V)| 2 ~ is observed. 
For noisy motion, (3=1, and the spectrum is contin- 
uous. For periodic motion (3 = 2 and the spectrum is 
discrete. For singular continuous spectrum (as in SNAs 
[23| | ) the scaling exponent satisfies 1 < (3 < 2 [23[. Here 
we find that the Fourier sum obeys scaling, with an expo- 
nent slightly greater than unity; see Figs.HD((a) and (b)). 
This implies that the dynamics in these attractors is typ- 
ically noisy unlike SNAs where the dynamical correlation 
persists over long times due to intermittency. For non- 
chaotic attractors in the Henon and the Lorenz systems, 
we find the exponents (3 = 1.05 and (3 = 1.04 respec- 
tively. Fig. EK(c) and (d)) shows the respective spectral 
trajectories in the complex plane (Re X,Im X). 

A detailed characterization of the nature of dynamics 
can be obtained from measures based on recurrences [24| . 
Recurrence plots (RPs) are defined for a given trajectory 
{*-i}iLi through the matrix 

Ri d = 0(5 x, ||), ij = 1, . . . ,N (23) 




Re X Re X 



FIG. 6: Finite-time Fourier power \X(Q, N)\ 2 versus T on a 
logarithmic scale for (a) the Henon map, and (b) the Lorenz 
system. Fractal walk of the spectral trajectory in the com- 
plex plane (ReX, ImX) for (c) the Henon map, and (d) for 
the Lorenz system. Parameter values are as in Fig. [Ha) and 
Fig. 121 



where S is a predefined threshold, O(-) the Heaviside 
function and || . || the maximum norm. The maximum 
norm (also called infinity norm) of a vector x of length 
N is given by || x ||oo— Tnax(\x\\, . . . , |xjv|)» Points that 
are closer (respectively further) than 5 yield an entry "1" 
(respectively "0") in the matrix Rij. Then, the values 
"1" and "0" are depicted as black and white dot in a 
two-dimensional plot, providing a visual representation 
of the system dynamics. The RPs exhibit characteris- 
tic large scale and small scale patterns (called typology 
and texture respectively); these have been comprehen- 
sively reviewed recently [25[. The selection criteria for 
the threshold 5 is discussed in details in a review by Mar- 
wan et al [25[. Here we take 6 in units of the standard 
deviation <j of the trajectory. 

In Fig. [7J we compare the RPs for a SNA in the 
quasiperiodically forced Henon map [27[ and the ANA 
for the modulated Henon map. It is clear that the RP 
of ANA (Fig. 0(b), (d)) consists of more isolated corre- 
lated points and short diagonal segments depicting short- 
range correlations. On the other hand, the RP for SNA 
(Fig. Ufa), (c)) has a larger distribution of longer diagonal 
line segments implying that correlation persists over long 
times - a signature of quasiperiodic driving. 
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FIG. 7: Comparision of recurrence plots for (a) a strange nonchaotic attractor for the quasiperiodically forced Henon map 
[27( ] . and, (b) the aperiodic nonchaotic attractor for the driven Henon map with parameters as given in Fig. [TJb). The RP for 
SNA has longer diagonal length segments due to long range correlations as compared to that of the ANA which consists of 
more isolated correlated points and shorter diagonal length segments. This is seen in (c) and (d) where the diagonal length 
distribution (with minimum length L m in = 2 in each case) for the two cases are compared. We calculated the determinism 
(DET) for both the cases (see Eq. 46 in Marwan et al 25]) and for SNA, BET = 0.78, whereas for ANA it is BET = 0.62. 
Here we take S = 0.3a where a is the standard deviation of the trajectory. 



EXPERIMENT 



In the electronic experiments reported here we have 
designed a circuit that essentially obeys the Lorenz equa- 
tions [17[ and permits one of the parameters to be 
switched between two values. 

The typical parameters for the butterfly attractor are 
<j = 10, 6 = 8/3 and in our circuit the value of p may be 
switched between 28 and 7 (as in Fig. [4|). The switch- 
ing of p is controlled by a chaotic pulse generated from 
a Chua circuit [28|, [29[: p = 28 if the drive signal has 
value 1, else p— 7 (and the drive signal has the value 
0). The experimental circuit of the pulse driven Lorenz 
oscillator is shown in Fig. [51 A Chua circuit is designed 
using two op- amps (U1-U2: /i741), capacitors CI and 
C2, inductor LI with a leakage resistance R8 and other 
resistances R1-R7. It generates a chaotic double scroll 
for choice components noted in the circuit diagram. The 
dynamics of the Chua circuit can be controlled by vary- 
ing Rl resistance keeping other components fixed. The 
double-scroll chaos from the Chua circuit is then applied 



to a Schmitt trigger circuit designed by using op-amp 
U3, an inverting amplifier U4 and associated resistances 
R9-R13. The output from U3 and U4 are used to con- 
trol the analog switches U5A and U6A respectively to 
allow continuity of either R16 or R17 in the Lorenz cir- 
cuit. The Lorenz circuit is implemented using two analog 
multipliers U7-U8, and three op-amps (U9-U11: /i741), 
capacitors C3-C5 and resistances R14-R21. The choice of 
resistances R16 and R17 made the selection of p- value be- 
tween 7 (= R19/R16) and 28 (= R19/R17) respectively. 
The other parameter of the Lorenz circuit are decided as 
a = R19/R14 and b = R19/R21. The analog switches 
are in ON state if their control pulse at VC terminal is 
positive. So the analog switch U5A is in ON state when 
the output of U3 is positive but U6A is in OFF state. 

Alternately, the analog switch U6A is in ON state and 
U5A is in OFF state when the output of U3 is negative 
but inverted by the U4 to make the control pulse positive 
at the VC terminal of U6A. The oscilloscope picture of 
the control pulse as generated from the Chua circuit is 
shown in Fig. [9l The upper trace is the double scroll 
chaotic signal from the Chua circuit which is processed 



7 




FIG. 8: (Colour Online) Experimental Lorenz circuit: ICs (U1-U8) power supply ±9Volt and (U9-U11) power supply ±15 
Volt. All resistances are in Ohms. 




FIG. 9: (Colour Online) Chaotic pulse generated using Chua 
oscillator: upper trace is the chaotic voltage VC2 at capacitor 
node C2 of the Chua circuit and the lower trace is the voltage 
measured at the output of the op- amp U3. Upper trace is 
scaled-up and the lower trace is scaled-down for visual clarity. 



by the Schmitt trigger U3. The chaotic pulse is clearly 
seen in the lower trace as switching between a positive 
and a negative value almost randomly; the signal is scaled 
down in the oscilloscope. 

The chaotic control signal switches the p-value of the 
Lorenz circuit aperiodically. The phase portrait of the 
Lorenz circuit is shown in Fig. [lOj this is the ANA that 
results from the theoretical strategy outlined in Sec HIT} 
see Fig. SJ 

To show that the attractors obtained experimentally 
satisfy the synchronization condition, we constructed an 




FIG. 10: (Colour Online) Oscilloscope picture of the ANAs in 
the modulated Lorenz circuit: output voltage of U8 plotted 
against the output voltage of Ull. 



auxiliary system [30[. Two identical Lorenz oscillators 
are controlled by the chaotic pulse generated from a sin- 
gle Chua oscillator. The components of the Lorenz os- 
cillators are carefully chosen with 1% tolerance so that 
both the oscillators are almost identical. Synchronization 
of two the nonchaotic Lorenz circuit is now investigated, 
the circuit scheme of which is shown in Fig. [TTJ We ob- 
serve that the output voltage of the Lorenz oscillators, 
OS-1 and OS-2, are completely synchronized. The oscil- 
loscope pictures of the two time series from the oscillators 
are shown in green and blue in Fig. [12] for comparison. 
The time series plotted one against the other, is shown as 
a thick line, confirms complete synchronization of the two 
oscillators (within experimental bounds). The width of 
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FIG. 13: (Colour Online) Oscilloscope picture of synchroniza- 
tion: similar output signals as shown in Fig. [TT] measured from 
two Lorenz circuits are plotted against each other to draw the 
synchronization manifold. 



FIG. 11: Block diagram of a system of two Lorenz circuits 
modulated by the chaotic drive of a Chua circuit. 




FIG. 12: (Colour Online) Oscilloscope picture of the out- 
put voltages showing the experimental time series of the two 
Lorenz circuits. 



the synchronization manifold is due to natural parameter 
mismatch between the two designed Lorenz oscillators; 
this is unavoidable in experiments. 



VI. DISCUSSION 

Motion that is both stable and aperiodic is ubiquitous 
in natural systems p], [3l| . The manner in which such dy- 
namics can be created is therefore of interest. One class 
of attractors that have these features has been known 
for some time now, but a quasiperiodic drive is essential 
for their creation |2j and thus these appear to be some- 
what exceptional. An area where these considerations 
are potentially important is in the dynamics of biological 
systems. Although not manifestly periodic, several bio- 



logical phenomena are stable, at least in a homeostatic 
sense [ail, [32|. Thus it is a moot question whether ape- 
riodic but nonchaotic attractors are responsible for such 
stability. 

On SNAs there is a delicate balance between global sta- 
bility: as was established by Sturman and Stark [33] there 
is an unstable set embedded within the attractor. The 
design strategy that we have enunciated in the present 
work keeps this feature in mind: the scheme we have pro- 
posed here is to modulate system parameters in such a 
manner as to achieve global stability while ensuring local 
instability. 

This method of dichotomous modulation creates at- 
tractors which are nonchaotic and have a fractal geom- 
etry on experimentally accessible timescales. We have 
recently shown [34| that SNAs created via quasiperiodic 
forcing are a manifestation of weak generalized synchro- 
nization, and that similar stable attractors can be created 
by chaotic forcing [34]. The parametric modulation used 
in the present case retains the skew-product structure 
of the dynamical system. The formation of these stable 
attractors can in some sense be seen as an instance of 
generalized synchronization [301 ]. 

The attractors created via such parameter modula- 
tion are quite distinct from SNAs. The Gaussian na- 
ture of the FTLE distribution shows that the dynamics 
is not intermittent. From the spectral properties it is 
evident that unlike SNAs, the power spectrum varies as 
\X(Q, N)\ 2 ~ TV, which occurs when the motion is ran- 
dom or chaotic. Our system being nonchaotic, the ran- 
domness in the motion comes from the stochastic nature 
of the modulation. This is confirmed by looking at the 
correlation properties via recurrence plots, which shows 
similiar behaviour to that of random or chaotic dynamics. 
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The aperiodic nonchaotic attractors can be realized in 
an experimental setup. As an example, we construct 
an electronic circuit experiment wherein a Chua double 
scroll attractor is used to drive a Lorenz attractor. The 
experiment closely matches the simulation result, and by 
experimentally constructing an auxiliary system [35[ we 
demonstrate that for a given drive sequence, trajectories 
with different initial states synchronize rapidly on the 
ANAs. This demonstrates the possibility of creating — or 
using — such dynamics in practical applications |36j . Fur- 



thermore such a realization shows the robustness of the 
proposed design scheme against external noise. 
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